Tumour Infiltrating Lymphocytes (TILs) and immune composition in breast cancer patients from Kenya: Spatial distributions and associations with risk factors and tumour characteristics.

BACKGROUND: The immune landscape of breast cancer (BC) in patients from Sub Saharan Africa is understudied. Our aims were to describe the distribution of Tumour Infiltrating Lymphocytes (TILs) within the intratumoural stroma (sTILs) and the leading/invasive edge stroma (LE-TILs), and to evaluate TILs across BC subtypes with established risk factors and clinical characteristics in Kenyan women. METHODS: Visual quantification of sTILs and LE-TILs were performed on Haematoxylin and eosin -stained pathologically confirmed BC cases based on the International TIL working group guidelines. Tissue Microarrays were constructed and stained with immunohistochemistry (IHC) for CD3, CD4, CD8, CD68, CD20, and FOXP3. Linear and logistic regression models were used to assess associations between risk factors and tumour features with IHC markers and total TILs, after adjusting for other covariates. RESULTS: A total of 226 invasive BC cases were included. Overall, LE-TIL (mean = 27.9, SD = 24.5) proportions were significantly higher than sTIL (mean = 13.5, SD = 15.8). Both sTILs and LE- TILs were predominantly composed of CD3, CD8, and CD68. We found higher TILs to be associated with high KI67/high grade and aggressive tumour subtypes, although these associations varied by TIL locations. Older age at menarche (≥ 15 vs. < 15 years) was associated with higher CD3 (OR: 2.06, 95%CI:1.26-3.37), but only for the intra-tumour stroma. CONCLUSION: The TIL enrichment in more aggressive BCs is similar to previously published data in other populations. The distinct associations of sTIL/LE-TIL measures with most examined factors highlight the importance of spatial TIL evaluations in future studies.

DNA methylation age in paired tumor and adjacent normal breast tissue in Chinese women with breast cancer.

BACKGROUND: Few studies have examined epigenetic age acceleration (AA), the difference between DNA methylation (DNAm) predicted age and chronological age, in relation to somatic genomic features in paired cancer and normal tissue, with less work done in non-European populations. In this study, we aimed to examine DNAm age and its associations with breast cancer risk factors, subtypes, somatic genomic profiles including mutation and copy number alterations and other aging markers in breast tissue of Chinese breast cancer (BC) patients from Hong Kong. METHODS: We performed genome-wide DNA methylation profiling of 196 tumor and 188 paired adjacent normal tissue collected from Chinese BC patients in Hong Kong (HKBC) using Illumina MethylationEPIC array. The DNAm age was calculated using Horvath's pan-tissue clock model. Somatic genomic features were based on data from RNA sequencing (RNASeq), whole-exome sequencing (WES), and whole-genome sequencing (WGS). Pearson's correlation (r), Kruskal-Wallis test, and regression models were used to estimate associations of DNAm AA with somatic features and breast cancer risk factors. RESULTS: DNAm age showed a stronger correlation with chronological age in normal (Pearson r = 0.78, P < 2.2e-16) than in tumor tissue (Pearson r = 0.31, P = 7.8e-06). Although overall DNAm age or AA did not vary significantly by tissue within the same individual, luminal A tumors exhibited increased DNAm AA (P = 0.004) while HER2-enriched/basal-like tumors exhibited markedly lower DNAm AA (P = < .0001) compared with paired normal tissue. Consistent with the subtype association, tumor DNAm AA was positively correlated with ESR1 (Pearson r = 0.39, P = 6.3e-06) and PGR (Pearson r = 0.36, P = 2.4e-05) gene expression. In line with this, we found that increasing DNAm AA was associated with higher body mass index (P = 0.039) and earlier age at menarche (P = 0.035), factors that are related to cumulative exposure to estrogen. In contrast, variables indicating extensive genomic instability, such as TP53 somatic mutations, high tumor mutation/copy number alteration burden, and homologous repair deficiency were associated with lower DNAm AA. CONCLUSIONS: Our findings provide additional insights into the complexity of breast tissue aging that is associated with the interaction of hormonal, genomic, and epigenetic mechanisms in an East Asian population.

Gene Expression Profiling Identifies Two Chordoma Subtypes Associated with Distinct Molecular Mechanisms and Clinical Outcomes.

PURPOSE: Chordoma is a rare bone tumor with a high recurrence rate and limited treatment options. The aim of this study was to identify molecular subtypes of chordoma that may improve clinical management. EXPERIMENTAL DESIGN: We conducted RNA sequencing in 48 tumors from patients with Chinese skull-base chordoma and identified two major molecular subtypes. We then replicated the classification using a NanoString panel in 48 patients with chordoma from North America. RESULTS: Tumors in one subtype were more likely to have somatic mutations and reduced expression in chromatin remodeling genes, such as PBRM1 and SETD2, whereas the other subtype was characterized by the upregulation of genes in epithelial-mesenchymal transition and Sonic Hedgehog pathways. IHC staining of top differentially expressed genes between the two subtypes in 312 patients with Chinese chordoma with long-term follow-up data showed that the expression of some markers such as PTCH1 was significantly associated with survival outcomes. CONCLUSIONS: Our findings may improve the understanding of subtype-specific tumorigenesis of chordoma and inform clinical prognostication and targeted options.

Rare germline variants in PALB2 and BRCA2 in familial and sporadic chordoma.

Chordoma is a rare bone tumor with genetic risk factors largely unknown. We conducted a whole-exome sequencing (WES) analysis of germline DNA from 19 familial chordoma cases in five pedigrees and 137 sporadic chordoma patients and identified 17 rare germline variants in PALB2 and BRCA2, whose products play essential roles in homologous recombination (HR) and tumor suppression. One PALB2 variant showed disease cosegregation in a family with four affected people or obligate gene carrier. Chordoma cases had a significantly increased burden of rare variants in both genes when compared to population-based controls. Four of the six PALB2 variants identified from chordoma patients modestly affected HR function and three of the 11 BRCA2 variants caused loss of function in experimental assays. These results, together with previous reports of abnormal morphology and Brachyury expression of the notochord in Palb2 knockout mouse embryos and genomic signatures associated with HR defect and HR gene mutations in advanced chordomas, suggest that germline mutations in PALB2 and BRCA2 may increase chordoma susceptibility. Our data shed light on the etiology of chordoma and support the previous finding that PARP-1 inhibitors may be a potential therapy for some chordoma patients.

Association of Genetic Ancestry With Terminal Duct Lobular Unit Involution Among Healthy Women.

Reduced age-related terminal duct lobular unit (TDLU) involution has been linked to increased breast cancer risk and triple-negative breast cancer. Associations of TDLU involution levels with race and ethnicity remain incompletely explored. Herein, we examined the association between genetic ancestry and TDLU involution in normal breast tissue donated by 2014 healthy women in the United States. Women of African ancestry were more likely than European women to have increased TDLU counts (odds ratio [OR]trend = 1.36, 95% confidence interval [CI] = 1.07 to 1.74), acini counts per TDLU (OR = 1.47, 95% CI = 1.06 to 2.03), and median TDLU span (ORtrend = 1.44, 95% CI = 1.08 to 1.91), indicating lower involution, whereas East Asian descendants were associated with decreased TDLU counts (ORtrend = 0.52, 95% CI = 0.35 to 0.78) after controlling for potential confounders. These associations are consistent with the racial variations in incidence rates of triple-negative breast cancer in the United States and suggest opportunities for future work examining whether TDLU involution may mediate the racial differences in subtype-specific breast cancer risk.

Integrated Analysis of Coexpression and Exome Sequencing to Prioritize Susceptibility Genes for Familial Cutaneous Melanoma.

The application of whole-exome sequencing has led to the identification of high- and moderate-risk variants that contribute to cutaneous melanoma susceptibility. However, confirming disease-causing variants remains challenging. We applied a gene coexpression network analysis to prioritize the candidate genes identified from whole-exome sequencing of 34 melanoma-prone families, with at least three affected members sequenced per family (N = 119 cases). A coexpression network was constructed from genotype-tissue expression project, skin melanoma from the cancer genome atlas, and primary melanocyte cultures. We performed module-specific enrichment and focused on modules associated with pigmentation processes because they are the best-studied and most well-known risk factors for melanoma susceptibility. We found that pigmentation-associated modules across the four expression datasets examined were enriched for well-known melanoma susceptibility genes plus genes associated with pigmentation. We also used network properties to prioritize genes within pigmentation modules as candidate susceptibility genes. Integrating information from coexpression network analysis and variant prioritization, we identified 36 genes (such as DCT, TPCN2, TRPM1, ATP10A, and EPHA5) as potential melanoma risk genes in the families. Our approach also allowed us to link families with private gene mutations on the basis of gene coexpression patterns and thereby may provide an innovative perspective in gene identification in high-risk families.

Comparison of somatic mutation landscapes in Chinese versus European breast cancer patients.

Recent genomic studies suggest that Asian breast cancer (BC) may have distinct somatic features; however, most comparisons of BC genomic features across populations did not account for differences in age, subtype, and sequencing methods. In this study, we analyzed whole-exome sequencing (WES) data to characterize somatic copy number alterations (SCNAs) and mutation profiles in 98 Hong Kong BC (HKBC) patients and compared with those from The Cancer Genome Atlas of European ancestry (TCGA-EA, N = 686), which had similar distributions of age at diagnosis and PAM50 subtypes as in HKBC. We developed a two-sample Poisson model to compare driver gene selection pressure, which reflects the effect sizes of cancer driver genes, while accounting for differences in sample size, sequencing platforms, depths, and mutation calling methods. We found that somatic mutation and SCNA profiles were overall very similar between HKBC and TCGA-EA. The selection pressure for small insertions and deletions (indels) in GATA3 (false discovery rate (FDR) corrected p < 0.01) and single-nucleotide variants (SNVs) in TP53 (nominal p = 0.02, FDR corrected p = 0.28) was lower in HKBC than in TCGA-EA. Among the 13 signatures of single-base substitutions (SBS) that are common in BC, we found a suggestively higher contribution of SBS18 and a lower contribution of SBS1 in HKBC than in TCGA-EA, while the two APOBEC-induced signatures showed similar prevalence. Our results suggest that the genomic landscape of BC was largely very similar between HKBC and TCGA-EA, despite suggestive differences in some driver genes and mutational signatures that warrant future investigations in large and diverse Asian populations.

Prognostic significance of tumor-infiltrating lymphocytes in premenopausal, luminal breast cancer treated with adjuvant endocrine therapy.

PURPOSE: Tumor-infiltrating lymphocytes (TILs) have strong prognostic value in triple-negative and HER2-enriched breast cancer, but their prognostic role in luminal breast cancer (LBC) is less clear. Here, we assessed the overall TIL levels and CD8+ T-cells in relation to the prognosis of LBC patients from China. METHODS AND RESULTS: A total of 596 patients with LBC who were premenopausal and treated with adjuvant endocrine therapy were included. Among them, 160 cases were evaluated for CD8 by immunohistochemical (IHC) staining. Whole-section hematoxylin and eosin and IHC staining were visually assessed for stromal TILs (sTILs), stromal CD8+ T-cells (sCD8), and intratumoral CD8+ T-cells (iCD8). Multivariable analyses were used to test the associations between TILs and disease-free survival (DFS) and overall survival (OS) with the adjustment for clinicopathologic characteristics and treatment. High sTILs (≥10%) were associated with high histologic grade (P<0.001), luminal B/HER2- (P<0.001), luminal B/HER2+ subtype (P=0.002), and high Ki67 expression (≥25%; P=0.014). Similar associations were observed for sCD8 but not for iCD8. While sTILs and sCD8 were not associated with either DFS or OS, the presence of iCD8 (≥1%) was associated with better DFS in both univariate (HR=0.51, 95% CI 0.26-0.96, P=0.042) and multivariate (HR=0.48, 95% CI 0.25-0.92, P=0.027) analyses. Similar but less significant associations were found for iCD8 and OS (adjusted HR=0.35, 95% CI 0.11-1.10, P=0.073). CONCLUSIONS: Among Chinese premenopausal patients with LBC, iCD8 demonstrated suggestive associations with favorable outcome. In contrast, although sTILs and sCD8 were associated with more aggressive tumor features, they did not appear to be associated with clinical outcome.

Quantitative Mammographic Density Measurements and Molecular Subtypes in Chinese Women With Breast Cancer.

Background: Studies investigating associations between mammographic density (MD) and breast cancer subtypes have generated mixed results. We previously showed that having extremely dense breasts was associated with the human epidermal growth factor receptor-2 (HER2)-enriched subtype in Chinese breast cancer patients. Methods: In this study, we reevaluated the MD-subtype association in 1549 Chinese breast cancer patients, using VolparaDensity software to obtain quantitative MD measures. All statistical tests were 2-sided. Results: Compared with women with luminal A tumors, women with luminal B/HER2- (odds ratio [OR] = 1.20, 95% confidence interval [CI] = 1.04 to 1.38; P = .01), luminal B/HER2+ (OR = 1.22, 95% CI = 1.03 to 1.46; P = .03), and HER2-enriched tumors (OR = 1.30, 95% CI = 1.06 to 1.59; P = .01) had higher fibroglandular dense volume. These associations were stronger in patients with smaller tumors (<2 cm). In contrast, the triple-negative subtype was associated with lower nondense volume (OR = 0.82, 95% CI = 0.68 to 0.99; P = .04), and the association was only seen among older women (age 50 years or older). Conclusion: Although biological mechanisms remain to be investigated, the associations for the HER2-enriched and luminal B subtypes with increasing MD may partially explain the higher prevalence of luminal B and HER2+ breast cancers previously reported in Asian women.

Rare Germline Variants in Chordoma-Related Genes and Chordoma Susceptibility.

BACKGROUND: Chordoma is a rare bone cancer with an unknown etiology. TBXT is the only chordoma susceptibility gene identified to date; germline single nucleotide variants and copy number variants in TBXT have been associated with chordoma susceptibility in familial and sporadic chordoma. However, the genetic susceptibility of chordoma remains largely unknown. In this study, we investigated rare germline genetic variants in genes involved in TBXT/chordoma-related signaling pathways and other biological processes in chordoma patients from North America and China. METHODS: We identified variants that were very rare in general population and internal control datasets and showed evidence for pathogenicity in 265 genes in a whole exome sequencing (WES) dataset of 138 chordoma patients of European ancestry and in a whole genome sequencing (WGS) dataset of 80 Chinese patients with skull base chordoma. RESULTS: Rare and likely pathogenic variants were identified in 32 of 138 European ancestry patients (23%), including genes that are part of notochord development, PI3K/AKT/mTOR, Sonic Hedgehog, SWI/SNF complex and mesoderm development pathways. Rare pathogenic variants in COL2A1, EXT1, PDK1, LRP2, TBXT and TSC2, among others, were also observed in Chinese patients. CONCLUSION: We identified several rare loss-of-function and predicted deleterious missense variants in germline DNA from patients with chordoma, which may influence chordoma predisposition and reflect a complex susceptibility, warranting further investigation in large studies.

Impact of breast cancer risk factors on clinically relevant prognostic biomarkers for primary breast cancer.

PURPOSE: In addition to impacting incidence, risk factors for breast cancer may also influence recurrence and survival from the disease. However, it is unclear how these factors affect combinatorial biomarkers for aiding treatment decision-making in breast cancer. METHODS: Patients were 8179 women with histologically confirmed invasive breast cancer, diagnosed and treated in a large cancer hospital in Beijing, China. Individual clinicopathological (tumor size, grade, lymph nodes) and immunohistochemical (IHC: ER, PR, HER2, KI67) markers were used to define clinically relevant combinatorial prognostic biomarkers, including the Nottingham Prognostic Index (NPI: combining size, grade, nodes) and IHC4 score (combining ER, PR, HER2, KI67). Odds ratios (ORs) and 95% confidence intervals (CIs) for associations between breast cancer risk factors and quartiles (Q1-Q4) of NPI and IHC4 were assessed in multivariable polytomous logistic regression models. RESULTS: Overall, increasing parity (ORtrend(95% CI) = 1.20(1.05-1.37);Ptrend = 0.007), overweight (OR(95% CI)vs normal = 1.60(1.29-1.98)), and obesity (OR(95% CI) vs normal = 2.12(1.43-3.14)) were associated with higher likelihood of developing tumors with high (Q4) versus low (Q1) NPI score. Conversely, increasing age (ORtrend(95% CI) = 0.75(0.66-0.84);Ptrend < 0.001) and positive family history of breast cancer (FHBC) (OR(95% CI) = 0.66(0.45-0.95)) were inversely associated with NPI. Only body mass index (BMI) was associated with IHC4, with overweight (OR(95% CI) vs normal = 0.82(0.66-1.02)) and obese (OR(95% CI) vs normal = 0.52(0.36-0.76)) women less likely to develop high IHC4 tumors. Notably, elevated BMI was associated with higher NPI irrespective of hormone receptor-expression status. CONCLUSIONS: Our findings indicate that factors affecting breast cancer incidence, particularly age, parity, FHBC, and BMI, may impact clinically relevant prognostic biomarkers with implications for surveillance, prognostication, and counseling.

Tumor-Associated Stromal Cellular Density as a Predictor of Recurrence and Mortality in Breast Cancer: Results from Ethnically Diverse Study Populations.

PURPOSE: Tumor-associated stroma is comprised of fibroblasts, tumor-infiltrating lymphocytes (TIL), macrophages, endothelial cells, and other cells that interactively influence tumor progression through inflammation and wound repair. Although gene-expression signatures reflecting wound repair predict breast cancer survival, it is unclear whether combined density of tumor-associated stromal cells, a morphologic proxy for inflammation and wound repair signatures on routine hematoxylin and eosin (H&E)-stained sections, is of prognostic relevance. METHODS: By applying machine learning to digitized H&E-stained sections for 2,084 breast cancer patients from China (n = 596; 24-55 years), Poland (n = 810; 31-75 years), and the United States (n = 678; 55-78 years), we characterized tumor-associated stromal cellular density (SCD) as the percentage of tumor-stroma that is occupied by nucleated cells. Hazard ratios (HR) and 95% confidence intervals (CI) for associations between SCD and clinical outcomes [recurrence (China) and mortality (Poland and the United States)] were estimated using Cox proportional hazard regression, adjusted for clinical variables. RESULTS: SCD was independently predictive of poor clinical outcomes in hormone receptor-positive (luminal) tumors from China [multivariable HR (95% CI)fourth(Q4) vs. first(Q1) quartile = 1.86 (1.06-3.26); P trend = 0.03], Poland [HR (95% CI)Q4 vs. Q1 = 1.80 (1.12-2.89); P trend = 0.01], and the United States [HR (95% CI)Q4 vs. Q1 = 2.42 (1.33-4.42); P trend = 0.002]. In general, SCD provided more prognostic information than most classic clinicopathologic factors, including grade, size, PR, HER2, IHC4, and TILs, predicting clinical outcomes irrespective of menopausal or lymph nodal status. SCD was not predictive of outcomes in hormone receptor-negative tumors. CONCLUSIONS: Our findings support the independent prognostic value of tumor-associated SCD among ethnically diverse luminal breast cancer patients. IMPACT: Assessment of tumor-associated SCD on standard H&E could help refine prognostic assessment and therapeutic decision making in luminal breast cancer.

Whole genome sequencing of skull-base chordoma reveals genomic alterations associated with recurrence and chordoma-specific survival.

Chordoma is a rare bone tumor with an unknown etiology and high recurrence rate. Here we conduct whole genome sequencing of 80 skull-base chordomas and identify PBRM1, a SWI/SNF (SWItch/Sucrose Non-Fermentable) complex subunit gene, as a significantly mutated driver gene. Genomic alterations in PBRM1 (12.5%) and homozygous deletions of the CDKN2A/2B locus are the most prevalent events. The combination of PBRM1 alterations and the chromosome 22q deletion, which involves another SWI/SNF gene (SMARCB1), shows strong associations with poor chordoma-specific survival (Hazard ratio [HR] = 10.55, 95% confidence interval [CI] = 2.81-39.64, p = 0.001) and recurrence-free survival (HR = 4.30, 95% CI = 2.34-7.91, p = 2.77 × 10-6). Despite the low mutation rate, extensive somatic copy number alterations frequently occur, most of which are clonal and showed highly concordant profiles between paired primary and recurrence/metastasis samples, indicating their importance in chordoma initiation. In this work, our findings provide important biological and clinical insights into skull-base chordoma.

Mental health care utilization and stigma in the military: comparison of Asian Americans to other racial groups.

Objective: To investigate race disparities in the US Military among Asian, White, Black, Native American and Other, seeking mental health care in the context of stigma defined by perceived damage to career.Design: Using 2008 survey data taken from US military personnel, mental disorders including depression, generalized anxiety disorder, suicidal ideation, suicidal attempt and post-traumatic stress disorder serious psychological distress (as defined in Kessler - 6), as well as seeking mental health care in past 12 months and stigma were dichotomized and weighted logistic regression models were used.Results: A significant race disparity existed in seeking mental health care when data were stratified by stigma and depression adjusted for demographic variables. Compared to Asians with depression that perceived stigma, Blacks were more likely to seek mental health care (OR with 95% confidence interval for Asians: 3.97[2.21, 7.15], Black: 9.25[6.02, 14.20], p < .005) adjusting for demographic variables. Similar results held for other mental disorders with the exception of suicide attempts and serious psychological distress. Compared to Asians with serious psychological distress who did not perceive stigma, only Whites were more likely to seek mental health care (OR for Asians: 3.27[2.15, 4.97], White: 6.47[4.60, 9.11], p < .005). Among those without a mental health disorder, regardless of the presence or absence of perceived stigma, there was no disparity between any two race groups in seeking mental health care.Conclusion: Among individuals having perceived stigma with mental health disorders, Asian American active-duty personnel may be less likely to use mental health care when compared to non-Asian peers.

Age-related DNA methylation in paired normal and tumour breast tissue in Chinese breast cancer patients.

Age-related DNA methylation is a potential mechanism contributing to breast cancer development. Studies of primarily Caucasian women have identified many CpG sites of age-related methylation in non-diseased breast tissue possibly driving cancer development over time. There is a paucity of studies involving Asian women whose ages at breast cancer onset are usually younger than Caucasians. We identified the 181 most consistent age-related methylation events in non-diseased breast tissue across published studies. Age-related methylation events were measured in adjacent normal and breast tumour tissue in an exclusively Asian population at the previously identified age-related methylation sites. Age-related methylation was found in 118 probes in adjacent normal breast tissue. Methylation of 99% of these sites was increased with age and predominantly located on CpG islands in promoter regions. To ascertain biological relevance to breast cancer, we focused on the 37 sites with overall higher methylation in tumour compared to adjacent normal samples. Some sites positively related to age, including AQP5 and CORO6, inversely correlated with gene expression. Several others have known involvement in suppression of carcinogenesis including GPC5 and SST, suggesting that perturbation of epigenetic regulation at these sites due to ageing may contribute to the progression of carcinogenesis. This study highlights an age-related methylation landscape in non-tumour tissue, consistent not just across studies, but also across different populations. We present candidate age-related methylation sites warranting further investigation as potential epigenetic drivers of breast cancer. They may serve as potential targets of site-specific demethylation intervention strategies for the prevention of age-related breast cancer.

Using whole-exome sequencing and protein interaction networks to prioritize candidate genes for germline cutaneous melanoma susceptibility.

Although next-generation sequencing has demonstrated great potential for novel gene discovery, confirming disease-causing genes after initial discovery remains challenging. Here, we applied a network analysis approach to prioritize candidate genes identified from whole-exome sequencing analysis of 98 cutaneous melanoma patients from 27 families. Using a network propagation method, we ranked candidate genes by their similarity to known disease genes in protein-protein interaction networks and identified gene clusters with functional connectivity. Using this approach, we identified several new candidate susceptibility genes that warrant future investigations such as NGLY1, IL1RN, FABP2, PRKDC, and PROSER2. The propagated network analysis also allowed us to link families that did not have common underlying genes but that carried variants in genes that interact on protein-protein interaction networks. In conclusion, our study provided an analysis perspective for gene prioritization in the context of genetic heterogeneity across families and prioritized top potential candidate susceptibility genes in our dataset.

Immune gene expression profiling reveals heterogeneity in luminal breast tumors.

BACKGROUND: Heterogeneity of immune gene expression patterns of luminal breast cancer (BC), which is clinically heterogeneous and overall considered as low immunogenic, has not been well studied especially in non-European populations. Here, we aimed at characterizing the immune gene expression profile of luminal BC in an Asian population and associating it with patient characteristics and tumor genomic features. METHODS: We performed immune gene expression profiling of tumor and adjacent normal tissue in 92 luminal BC patients from Hong Kong using RNA-sequencing data and used unsupervised consensus clustering to stratify tumors. We then used luminal patients from The Cancer Genome Atlas (TCGA, N = 564) and a Korean breast cancer study (KBC, N = 112) as replication datasets. RESULTS: Based on the expression of 130 immune-related genes, luminal tumors were stratified into three distinct immune subtypes. Tumors in one subtype showed higher level of tumor-infiltrating lymphocytes (TILs), characterized by T cell gene activation, higher expression of immune checkpoint genes, higher nonsynonymous mutation burden, and higher APOBEC-signature mutations, compared with other luminal tumors. The high-TIL subtype was also associated with lower ESR1/ESR2 expression ratio and increasing body mass index. The comparison of the immune profile in tumor and matched normal tissue suggested a tumor-derived activation of specific immune responses, which was only seen in high-TIL patients. Tumors in a second subtype were characterized by increased expression of interferon-stimulated genes and enrichment for TP53 somatic mutations. The presence of three immune subtypes within luminal BC was replicated in TCGA and KBC, although the pattern was more similar in Asian populations. The germline APOBEC3B deletion polymorphism, which is prevalent in East Asian populations and was previously linked to immune activation, was not associated with immune subtypes in our study. This result does not support the hypothesis that the germline APOBEC3B deletion polymorphism is the driving force for immune activation in breast tumors in Asian populations. CONCLUSION: Our findings suggest that immune gene expression and associated genomic features could be useful to further stratify luminal BC beyond the current luminal A/B classification and a subset of luminal BC patients may benefit from checkpoint immunotherapy, at least in Asian populations.

Clinicopathological and epidemiological significance of breast cancer subtype reclassification based on p53 immunohistochemical expression.

TP53 mutations are common in breast cancer and are typically associated with more aggressive tumor characteristics, but little is known about the clinicopathological and epidemiological relevance of p53 protein expression, a TP53 mutation surrogate, in breast cancer subtypes. In this study of 7226 Chinese women with invasive breast cancer, we defined breast cancer subtypes using immunohistochemical (IHC) measures of hormone receptors and HER2 in conjunction with histologic grade. p53 expression status was then used to further stratify subtypes into p53-positive and p53-negative. Odds ratios (ORs) and 95% confidence intervals (CIs) in case-only logistic regression analyses were used to examine heterogeneity across different subtypes. The frequency of p53 protein expression varied by breast cancer subtype, being lowest in the luminal A-like and highest in the triple-negative and HER2-enriched subtypes (P-value < 0.01). In luminal A-like and B-like/HER2-negative subtypes, p53 positivity was associated with early-onset tumors, high grade, high proliferative index, and basal marker (CK5/6 and EGFR) expression. Further, compared with luminal A-like/p53-negative patients, A-like/p53-positive patients were more likely to be parous [adjusted OR parous vs. nulliparous = 2.67 (1.60, 4.51); P-value < 0.01] and to have breastfed [adjusted OR ever vs. never = 1.38 (1.03, 1.85); P-value = 0.03]. p53 positivity was not associated with examined clinical and risk factors in other tumor subtypes. Overall, these findings suggest that p53 expression, which is readily available in many settings, can be used to identify phenotypes of luminal A-like breast cancer with distinct clinical and epidemiological implications.

Associations between mammographic density and tumor characteristics in Chinese women with breast cancer.

PURPOSE: Mammographic density (MD) is a strong risk factor for breast cancer, yet its relationship with tumor characteristics is not well established, particularly in Asian populations. METHODS: MD was assessed from a total of 2001 Chinese breast cancer patients using Breast Imaging Reporting and Data System (BI-RADS) categories. Molecular subtypes were defined using immunohistochemical status on ER, PR, HER2, and Ki-67, as well as tumor grade. Multinomial logistic regression was used to test associations between MD and molecular subtype (luminal A = reference) adjusting for age, body mass index (BMI), menopausal status, parity, and nodal status. RESULTS: The mean age at diagnosis was 51.7 years (SD = 10.7) and the average BMI was 24.7 kg/m2 (SD = 3.8). The distribution of BI-RADS categories was 7.4% A = almost entirely fat, 24.2% B = scattered fibroglandular dense, 49.4% C = heterogeneously dense, and 19.0% D = extremely dense. Compared to women with BI-RADS = A/B, women with BI-RADS = D were more likely to have HER2-enriched tumors (OR = 1.81, 95% CI 1.08-3.06, p = 0.03), regardless of menopausal status. The association was only observed in women with normal (< 25 kg/m2) BMI (OR = 2.43, 95% CI 1.24-4.76, p < 0.01), but not among overweight/obese women (OR: 0.98, 95% CI 0.38-2.52, p = 0.96). CONCLUSIONS: Among Chinese women with normal BMI, higher breast density was associated with HER2-enriched tumors. The results may partially explain the higher proportion of HER2+ tumors previously reported in Asian women.